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Abstract. System analysis of oxygen production by hydrogen reduction of lunar regolith has shown the importance of the 
relative time scales for regolith heating and chemical reaction to overall performance. These values detennine the sizing 
and power requirements of the system and also impact the number and operational phasing of reaction chambers. In this 
paper, a Nusselt number correlation analysis is perfonned to detennine the heat transfer rates and regolith heat up times 
in a fluidized bed reactor heated by a central heating element (e.g., a resistively heated rod, or a solar concentrator heat 
pipe). A coupled chemical and transport model has also been developed for the chemical reduction of regolith by a 
continuous flow of hydrogen. The regolith conversion occurs on the surfaces of and within the regolith particles. Several 
important quantities are identified as a result of the above analyses. Reactor scale parameters include the void fraction 
(i.e., the fraction of the reactor volume not occupied by the regolith particles) and the residence time of hydrogen in the 
reactor. Particle scale quantities include the particle Reynolds number, the Archimedes number, and the time needed for 
hydrogen to diffuse into the pores of the regolith particles. The analysis is used to determine the heat up and reaction 
times and its application to NASA’s oxygen production system modeling tool is noted. 
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INTRODUCTION 

Hydrogen reduction of lunar regolith has been proposed as a viable technology for oxygen production on the moon. 
(Taylor and Carrier 1992, McKay and Allen 1996, Rice, Hermes and Musbah 1997). The reduction processes 
involve solid-gas reactions. Such systems are well developed for earth applications but mainly for the purpose of 
metal extraction as opposed to water vapor and oxygen production for lunar applications. In addition, in the lunar 
context, the reactor must interface with the other sub-system processes such as upstream regolith extraction and 
beneficiation, and downstream electrolysis and phase separation in a way that establishes the most favorable balance 
between efficiency, robustness, and equivalent system mass. For example, the reduction of ilmenite may be 
represented as: 


FeO.Ti0 2 + H 2 <=> Fe + Ti0 2 + H 2 0 . ( 1 ) 

It is carried out in the reactor sub-system. The formed water may be electrolyzed in the electrolyzer subsystem to 
produce oxygen and hydrogen. 

From upmass and cost considerations, it is essential to operate these systems on the lunar surface in as nearly a 
closed-loop manner as possible for the hydrogen. It is in this context that knowledge of processing times in the 
various sub systems becomes important. For example, a mismatch between the flux of water exiting the reactor and 
the amount of water that may be processed by the electrolyzer must be mitigated by design such as sizing the 
appropriate buffer tank. This paper investigates the factors that influence the overall processing time in the reactor 
sub system. 

The temperature of the reaction is one of the important parameters of the process. Ilmenite is just one constituent of 
mare basalt and it cannot be assumed that regolith beneficiation (e.g., in this case, removal of the non-ilmenite 
components) will be carried out prior to the reactor processes. The sintering and melting temperature of the 
unbeneficiated regolith is noted to be in the range of 1370-1670 K (Chambers 1995, McKay and Allen 1996). 
However, hydrogen reduction has been found to be effective also at lower temperatures, e.g., in the range of 970- 
1270 K (Yoshida et al., 2000). Therefore, while three-phase (gas-liquid-solid) reactors for hydrogen reduction may 



be applicable at post-melt temperatures, these are not considered here and attention is restricted to gas-solid 
interactions. In order to provide a basis for discussion, it is assumed that the chemical reduction of the regolith 
occurs in a gas-solid fluidized bed reactor. 

The conversion process in the reactor consists of the following overall steps: (a) transfer of regolith into the reactor, 
(b) fluidization and heating of the regolith, (c) reaction of hydrogen with the ilmenite portion of the regolith, (d) 
outflow of left-over hydrogen and produced water vapor, and (e) removal of the spent regolith from the reactor. In 
this paper, attention is focused on steps (b) -(d). 

The reported values for the thermal diffusivity, (1. of lunar regolith are very low (Heiken, Vaniman and French 
1991). In the absence of a gas flow, the time to heat the regolith bed over a distance 5 is proportional the ratio 5 2 /p, 
which results in unacceptably large values. On the other hand, hydrogen has a relatively high thermal conductivity 
and helps to stir the bed so that having hydrogen flow through the bed decreases the heating time appreciably. Thus, 
it is advantageous to fluidize the bed during the regolith heating step. There are also advantages related to reducing 
regolith particle agglomeration and sintering which can adversely impact reaction times and water production rates. 

The objectives of this paper are (a) to develop analytical models for regolith heating and for chemical reduction by 
hydrogen in order to identify parameters influencing the heating and reaction times, (b) to identify the parameters 
that govern these times, and (c) to obtain relationships for the reaction time in the two limits of diffusion (i.e., mass 
transfer) control and chemical equilibrium control. As noted earlier, these times are important from the overall 
oxygen production system perspective. The thermal analysis is based upon a Nusselt number correlation. The 
reaction model considered is based upon a “shrinking core” diffusion-reaction model for the gas-solid reduction 
process (Levenspiel 1998) which is coupled with the inflow/outflow of gases. The methodology has been discussed 
in our earlier work (Hegde, Balasubramaniam and Gokoglu 2007). It appears consistent with results of regolith 
particulate analysis conducted on Apollo-returned samples (McKay and Allen 1996) which suggest that the ilmenite 
particles are not completely solid but have some porosity that would enable gases to diffuse into the particles, i.e., 
the FeO portion can be essentially completely converted leaving behind the reduced iron. 

In the next section, the regolith heating is considered first and a model for the heating time and required power is 
developed. Pertinent results from the shrinking core analysis are also presented. Subsequently, in the Results 
section, factors impacting these processes are evaluated and discussed and applications of the model to NASA’s 
oxygen production system modeling tool are noted. 


MODEL DEVELOPMENT 

In the first sub-section below, the time taken to heat the regolith is calculated using a Nusselt number correlation for 
fluidized bed heat transfer. This is followed by a summary of the particle conversion analysis and the coupling 
between the particle conversion and the bulk reactor flow. 

Heat Transfer Analysis 

It is assumed that the regolith is heated by a cylindrical rod located centrally in the cylindrical reactor. The surface 
temperature of the heater rod can be controlled e.g., by adjusting the current through it as in an electrical heater or 
by controlling the energy flux as in a solar concentrator heat pipe. The different pathways for heat transfer in the bed 
are gas convection and conduction, particle motions, and possibly, radiative transfer. 

Fluidized bed heat transfer correlations developed in the literature under convective-conductive conditions suggest 
that the controlling step is the heat transfer from the heater surface to the gas film next to it. Direct transfer of heat to 
the bed particles adjoining the heater is not efficient because the contact area between particles and the heater 
surface is small. If the bed is well stirred by fluidization, its temperature is nearly uniform. In such a case, a Nusselt 
number may be defined which relates the heat transfer to the bed through the gas film in contact with the heater, i.e., 
the regolith bed is taken to be a uniform entity. The reactor wall is not usually insulated from the bed, so in general, 
there will also be heat transfer from the bed to the reactor wall. If the reactor wall is not too thick and its thermal 
conductivity is high, the wall may be assumed to be at a uniform temperature. Furthermore, if the reactor’s thermal 
mass (i.e., the product of its mass and specific heat) is small compared to that of the bed, it is reasonable to assume 
that the temperature of the reactor is the same as that of the bed. 



The heating rate with the above scenario is given by: 


M * r C dr ^ = hA dr (T l -T), 


dt 


(2) 


where the thin film heat transfer coefficient, h, is related to the Nusselt number, Nu, by: 

Nil = 2hr p / X . (3) 

In Eq.(2), M eff C e ff is the effective value of the product of the mass and specific heat and is a weighted average of 
the values for the gas, regolith, and reactor. A e ff is the effective or active surface area of the heater rod in contact 
with the bed. 

The solution of Eq. (2) for the instantaneous bed temperature T is given by: 

T = T h~ ( T i, ~ T i ) exp (~hA eff t !(M eff C eff )) . (4) 


Conversely, the time to raise the bed temperature from Tj, to a temperature T 2 is given by: 
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Various fluidized bed heat transfer Nusselt number correlations are available in the literature. For example, the 
following one (Borodulya et. al. 1991) appears to describe the situation adequately: 

Nu = 0.74^r 01 (^) 0 - 14 (-^)°- 24 (l-^) 23 + 0.46 Re Pr — — . (6) 

P s C p S £ 

The Archimedes number, Ar, the Reynolds number. Re, and the Prandtl number, Pr, are defined as follows: 

Ar = 8 gr]p g (p s - p g ) / // 2 , (7) 


Re = 2ur p p g / p , (8) 

Pr = C pg p /A g . (9) 

The energy, E, required to heat the bed and reactor from Tj to T 2 , in the absence of other losses (e.g., radiative), is 
given by: 


E = M eff C eff {T 2 -T x ). (10) 

An average power equired is obtained by dividing E by the time, t h . Note that the instantaneous power, 
MefAffDT/dt, can be obtained by substituting the solution for T (Eq. (4 )) into Eq. (2). 

Since radiative losses vary as the fourth power of temperature, they can become important as the reactor is heated to 
higher temperatures. It is expected that there would be insulation around the reactor that would reduce this loss. This 
loss can be accounted for by including it in the energy required. 



Chemical conversion analysis 


The chemical conversion analysis couples the “shrinking core” model with the inflow and outflow of the gases and 
has been discussed in earlier work (Hegde, Balasubramaniam and Gokoglu, 2007). Here, some results are presented 
to highlight the findings from the model and to show how the characteristic times of the process are obtained. 

Defining the conversion a by: 

r 3 

a = l-^r, (ll) 

r P 

the shrinking core analysis of a single particle yields a relationship between a and the associated conversion time, t: 


3[1 -«] 2/3 -2[l-a]-l = - 
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( 12 ) 


Further, assuming equilibrium conditions at the reaction surface, the above equation reduces to: 

6 DkcJ 


3[l-a] 2/3 - 2[l-a]-l = - 


(1 + k) Pa r 2 p 


(13) 


Note that an expression for the equilibrium constant, k, as a function of only the reaction temperature has been 
provided (William and Erstfeld, 1979): 


_ i n (-2126.1/7’+0.6439) 


k = 10 


(14) 


The time for complete conversion is given by setting a to unity. 

Coupling Particle Analysis with the Reactor Bulk Flow 

Assuming that the gas velocity, concentration of hydrogen and the associated conversion a are functions of axial 
distance along the reactor, and the inlet hydrogen concentration is constant in time, the following equation has been 
obtained: 
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where a 0 is the value of a at x=0 and is given by Eq. (13). 

In certain cases, closed form solutions may be obtained for the concentration of hydrogen and extent of conversion 
in the reactor. An example is when the regolith and gases are sufficiently stirred or well-mixed so that gradients in 
conversion and gas concentration are smeared out. In this well-stirred limit, the concentration of hydrogen exiting 
the reactor is also the average concentration in the reactor and is given by: 
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RESULTS 


Regolith batch sizing will depend upon the oxygen production required and the processing time per batch. In 
general, the processing time per batch is the sum of the heating and chemical conversion times. This section 
considers factors influencing these times based upon the analysis presented earlier. 

Heating Time 

Equation (5) indicates that the heating time is related directly to the Nusselt number and batch size (i.e., regolith 
mass) and inversely related to the active surface area of the heater. The Nusselt number depends upon the Reynolds 
number of the gas flow relative to the regolith particles, the Archimedes number, and the Prandtl number of the 
fluidizing gas. The Prandtl number is practically constant with temperature. The Reynolds number is related to the 
fluidizing velocity and particle size. The particle size is fixed by the regolith characteristics and there may be some 
flexibility in this regard depending upon the preprocessing facilities available. The fluidizing velocity is also related 
to the particle size. The Archimedes number is mainly dependent upon the regolith density (which is much greater 
than the gas density). The most flexible way of changing the heating time seems to be by changing the regolith batch 
size, active heater area, and/or power input to the heater. 

The active heater surface area may be increased by increasing the heater diameter. For a fixed reactor diameter, this 
also increases the height, L, of the regolith bed which is also conducive to heat transfer. However, this may have 
implications for the reaction time as will be discussed later. The regolith batch size may be reduced to decrease the 
heating time but this may be restricted from overall oxygen production considerations. 

The power input to the heater is related to the heater temperature T h . From Eq. (5), it is clear that increasing the 
heater temperature will decrease the time to heat the regolith to a given temperature T 2 . However, there are 
limitations to raising Th substantially over T 2 e.g., the presence of hot spots or regions in the regolith bed which are 
at higher than desired temperature. 


Chemical Time 

The shrinking core analysis provides an expression for the conversion time of a single particle. In many instances it 
is assumed that the conversion time of the regolith bed is the same as the single particle conversion time. Equation 
(15) provides a means for determining the conditions under which this holds. Choosing x = L (i.e., the height of the 
regolith bed), it is clear that the right hand side of Eq. (15) is related to the ratio of the gas residence time in the bed, 
L/u, to its diffusion time in a particle, r p 2 /D. This time factor ratio plays an important role in determining the extent 
of conversion in the bed. For example, the upstream end (entrance region) of the bed is always exposed to fresh 
hydrogen and would be the first to get converted. The time ratio factor gives a measure of the conversion in the 
remaining portions of the bed. Figure 1 plots the conversion, a(L), at the downstream end (exit region) as a function 
of the time ratio factor for the case where do is equal to unity, i.e., when the particles at the entrance are already 
completely converted. As the time ratio factor becomes small (i.e., the intra-particle diffusion time is much greater 
than the gas residence time in the regolith bed) it can be seen from Fig. 1 that the conversion, a(L), also approaches 
unity, i.e., the particles at the downstream end of the bed (and more generally, at any location in the bed) are also 
consumed fully. This indicates that the time rate of chemical conversion for all particles in the bed is given by Eq. 
(13 ). The corresponding time for complete conversion may be termed the chemical time in the diffusion controlled 
limit. Setting a to one in Eq. (13) yields: 
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FIGURE 1 . Dependence of the Regolith Conversion at the Downstream End of the Reactor Bed Upon the Time 
Ratio Factor 3(l-s) DL/r p 2 u. 

On the other hand, as the residence time becomes much larger compared to the diffusion time, i.e., the time ratio 
factor becomes large, it can be seen from Fig. 1 that the conversion occurs first primarily near the upstream end of 
the bed and not much conversion occurs downstream since a(L) is small for this condition. This is referred to as the 
chemical equilibrium limit, i.e., the residence time in the bed is long enough so that the water concentration reaches 
its chemical equilibrium value in the gas phase. At this point no further conversion takes place. In this limiting case, 
the chemical time may be estimated by assuming equilibrium concentration of water in the outflow. Then, the molar 
outflow of water is given by : 


and the total molar water yield is: 
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(19) 


Equating this to the available molar water yield in the regolith provides the chemical time in the residence time limit 
or equilibrium limit: 


tResUmi, = Pad ~ £ ) V bed ' ¥ ■ ( 20 ) 

Since the conversion time in the diffusion limit depends inherently upon the particle characteristics such as its size 
and the intra-particle gas diffusivity, it is clear that this is a lower limit to the conversion time. The equilibrium 


limiting conversion time is variable however, depending upon the flow velocity and bed height. Varying these 
quantities is a way to slow down the conversion rate. This may be useful from an overall system design perspective. 

Application to System Design 

The analysis of heating and chemical times presented here may be applied to sizing and design of the reactor 
assembly. It has been noted in earlier studies that the regolith reduction rate increases directly with the temperature, 
(McKay and Allen 1996) presumably due to increase in the effective diffusion coefficient (Hegde, Balasubramaniam 
and Gokoglu 2007). This implies that most of the conversion occurs at the operating temperature. In other words, the 
heating and chemical reduction are essentially serial operations for a given batch of regolith. Thus, a system 
comprising two reactors may be used for increased throughput. Ignoring for current purposes, the regolith loading 
and unloading times into the reactor chamber, it is clear that designing for equal heating and reaction times will 
ensure continuous production of oxygen, i.e., as one batch of regolith is being heated the other batch is being 
chemically reduced. 

The equations provided here provide a framework for achieving the proper design parameters. The diffusion limited 
chemical time (Eq. (17)) provides the minimum time for conversion. Equations (5) and (10) in conjunction provide 
the power and heater rod characteristics (e.g., effective area) for a given heating time, e.g., the minimum reaction 
time. If the power available is less than the heating power required, the heating time has to be increased. In such a 
case, the reaction time needs to be modified e.g., by driving the reactor dimensions toward the residence limited 
conversion time (Eq. (20)). Clearly, this is an iterative process. Currently, NASA is developing an overall system 
model incorporating these ideas (Santiago-Maldonado 2007). 


CONCLUSIONS 

Analytical models for regolith heating and for chemical reduction by hydrogen have been developed to identify 
parameters influencing the heating and reaction times. These times are important from the overall oxygen 
production system perspective. The regolith heating model is based upon a Nusselt number correlation for fluidized 
bed heat transfer. While particle size and parameters such as the particle Reynolds and Archimedes numbers play a 
role in the heating time determination, it appears that the most feasible way of controlling this time is by adjusting 
the active heater surface area and mass of the regolith per batch. The chemical reduction model results from 
coupling a shrinking core type of analysis with the overall regolith bed flow characteristics. It is shown that the 
chemical time in the diffusion limit is also the lower limit for the conversion time. However, by adjusting the 
residence time of the gases in the regolith bed, the conversion time can be increased. 


NOMENCLATURE 


a = particle conversion 

a Q = particle conversion at x =0 

P = regolith thermal diffusivity (m 2 /s) 

5 = thermal conduction length scale (m) 

e = void fraction in regolith bed 

p a = molar concentration of ilmenite in regolith 

(mol/m^) 

Pg = gas density (kg/m^) 

p s = regolith density (kg/m^) 

p = gas viscosity (Ns/m~) 

X = gas thermal conductivity 
T = molar outflow of water (mol/s) 

A = reactor cross section area (m 2 ) 

A e ff = effective surface area of heater (m 2 ) 


Ar = Archimedes number 

c = molar concentration of hydrogen in particle 
(moles/m 3 ) 

c s = molar concentration of hydrogen at reaction 
surface (moles/m^) 

c x = molar concentration of hydrogen at location 
x in reactor (moles/ nr"’) 

c 0 = molar concentration of hydrogen at upstream 

end of regolith bed (moles/nr*) 

C e fi= effective specific heat (J/kg/K) 

Cpg = specific heat of gas (J/kg/K) 

Cps = specific heat of regolith (J/kg/K) 

D = effective gas diffusion coefficient (m-/s) 

E = energy required to heat bed 
h = heat transfer coefficient 



g = gravitational acceleration (m/s^) 

k = equilibrium constant 

L = length of regolith bed (m) 

M eff = effective mass (kg) 

M tota i = total molar water yield (mol) 

N = number of regolith particles in the bed 
Nu = Nusselt number 

Pr = Prandtl number 

i p = regolith particle radius (m) 

Re = Reynolds number 

t = time (s) 


t c = chemical time (s) 

tpj = intra-particle diffusion time (s) 

t D iimit = chemical time in diffusion limit (s) 
tR = reactor residence time (s) 

t R e S Limit = chemical time in equilibrium limit (s) 

T = temperature (K) 

T i, = heater temperature (K) 

Tj, T 2 = initial and final temperature of bed (K) 
u = velocity through reactor (m/s) 

Vbed = volume of regolith bed (m 3 ) 

x = axial coordinate (m) 
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